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We present one inherent shortcoming of the LDA+U method in respect of its self-interaction 
correction of the LDA. By reexamining the mean-field approximation on the Hubbard energy in 
the Hartree-Fock form, we have derived a new expression for the "double-counting" energy which 
led to a more reasonable self-interaction correction in an alternative LDA+U scheme. For the bulk 
Gd metal, the new scheme resulted in electronic properties more consistent with experiments in 
comparison to the LDA and the existing LDA+U schemes. 



Since its advent more than two decades ago [TJ [2], 
the LDA+U method has remained as one of the most 
efficient methods that describe the strong correlations 
of electrons confined in an atomic region. (Through- 
out this paper LDA includes its generalization to spin- 
polarized systems.) When applied to the transition metal 
and rare-earth metal compounds, the LDA+U method 
gives a significant improvement compared with the LDA 
not only for excited-state properties but also for ground- 
state properties. However, further applications of the 
method may require its improvement. For example, it 
has been pointed out that weakly-correlated systems are 
not properly described by the LDA+U and that a cor- 
rect prescription of the "double-counting" energy may 
be essential |5]. For strongly-correlated systems, there 
can exist drastic discrepancies in the locations of excited 
Hubbard states in comparison to the self-interaction- 
corrected LDA (SIC-LDA) calculations |4], which war- 
rants a closer look at the self-interaction correction incor- 
porated in the LDA+U. Herein, we present a shortcom- 
ing associated with its self-intereaction correction to the 
LDA, and provide an alternative LDA+U scheme with 
a new double-counting energy expression that does not 
bear the same problem, producing the results consistent 
with the SIC-LDA calculations. 

The LDA+U total energy functional takes the form 

E total [p a ,h] = E LDA [ Pa ]+E u (h) , (1) 

where Elda\Po\ is the usual local spin-density functional 
of the total electron spin densities (p CT (r), a =f, J,). The 
correction term to the LDA energy is given as 

E u {n) = E^(n)-Et DA (h) . (2) 

E ee is an electron-electron interaction energy of corre- 
lated d(f) electrons in orbitals {4> m }, typically described 
by the multiband Hubbard model in the (unrestricted) 
Hartree-Fock (HF) approximation with a density matrix 
n El- EldA ls t ne double-counting term that accounts 
approximately for the LDA-type energy of the correlated 
electrons. When electrons are in {cf> m }, they feel the extra 
potential, v J ^ na , from E u (h) in addition to the LDA po- 
tential. The effective potential for those correlated elec- 



trons is given as 

"m/( r ) = v ne(r) + v H (r) + v xca (r) + v^ na . (3) 
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FIG. 1: ls spin-orbital energy changes of a hydrogen atom by 
subsequent inclusions of vh = U, {v xc f ~ — J, v xc i ~ 0} and 
Vim — -f^(U—J) in the AL scheme. The dashed levels are the 
unoccupied j-spin energy level after the exchange splitting. 

Depending on the way the Ef^ DA expression is de- 
rived, the current LDA+U methods are classified into two 
schemes, AMF (around mean field) [T] and AL (atomic 
limit) (2]. A drawback of the AMF scheme has been rec- 
ognized, as it provides no effect on half-filled fully spin- 
polarized systems. The AL scheme does not have such 
a problem and hence is a preferred method for half-filled 
electron systems. However, we may recognize its short- 
coming when we consider its v J ^ la potential, 

v ^{r) = -{U-J){n mlT -^j . (4) 

An illustrating case is an isolated hydrogen atom, a half- 
filled system. With N = 1, ni s j = 1 and ni s ^ = 0, 
vf//(r) of Eq. [3] is given as 

= "ne(r) + U + v xca (r) + , (5) 

where the + and — signs are for the ls \- and |-spin 
orbitals of the hydrogen atom. Figure 1 shows how the 
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Is spin-orbitals change their energy, as vh = U, {v xc ^ w 
— J, v xc i w 0} and u* 1 ^. = + |(£/ — ^) are subsequently 
added. While the v^ a clearly corrects the LDA term, 
v xc ^ « —J, it does not bring the "f-spin orbital back 
to the original position. In practice, the (U — J) term 
in vf-^ may be parametrized into U' . However, the U' 
value that completes the correction of vh in Fig. 1 will 
provide an erroneously large exchange splitting energy 
(on-site repulsion energy, in this case). This problem can 
be removed in an alternative LDA+U scheme proposed 
in the following. 

We start with the electron-electron interaction energy 
within the HF approximation by using a diagonalized 
density matrix n [l][2]: 

1 V" 

= ^ / j ^ mm 1 ^ina^in 1 '—a 
m,m f ,<r 

2 ^ ^ (Um,7n' "^mm / )^mcr'^m'cr • (6) 

Eq. [6] is reexpressed to collect all Coulomb interaction 
terms in the Hartree energy component Ejj at the begin- 
ning: 
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The t7 and J parameters are by definition averaged 
"self-orbital" Coulomb (or exchange) integral and "pair- 
orbital" exchange integral, respectively, for which a ro- 
tational invariance is assumed. The Eh term is often 
approximated as J U N 2 under the rotational invariance, 
where U is the usual average Coulomb integral. For a 
more rigorous treatment, however, we leave Eh as it is 
without the approximation. Alternatively, averaging of 
all the Coulomb (exchange) integrals gives the U (J) pa- 
rameter in the conventional approximate expression for 
the Ehf in Eq.|6] however, not with \UN 2 but with Eh'- 
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(8) 



Following the mean-field (MF) treatment in the AMF 
scheme [T] [3] , we obtain the LDA-type energy corre- 



sponding to Eq. [7] upon replacing all n mcr 's by equal oc- 
cupation numbers, jfp[, among the (21 + l)-degenerate 
orbitals: 
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The MF treatment may also affect numerical values of 
the U and J parameters, and hence the parameters are 
designated in Eq. [j)]as U MF and J MF , respectively. From 
Eq. [9] the spin-polarization energy can be estimated for 
an open-shell system with an M number of unpaired elec- 
trons by noting J2 a N * = \(N 2 + M 2 ): 
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(10) 



The quadratic dependence of the spin-polarization energy 
on the spin magnetic moment M is consistent with the 
previous studies within LDA j6]. In particular, our recent 
work has shown quantitatively |7] that within a second- 
order perturbational theory, the corresponding LDA ex- 
pression is given for isolated atoms as 
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(11) 



in which the U xc and J xc are defined as self-orbital and 
pair-orbital exchange-correlation interaction energies, re- 
spectively. They define the atomic exchange parame- 
ter [3, 



U xc + 2U X 
2/+ 1 



(12) 



which is analogous to the Stoner exchange parameter of 
bulk metals. The two different sets of exchange param- 
eters, calculated in LDA, are numerically identical for 
the 3d elements |7], which is expected to hold for the 4/ 
elements as well. 



From the equivalency of the expressions in Eqs. 10 and 
[TT| we find the following relationship in replacement of 
Eq.[9] 
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(13) 



The second term in Eq. 13 can be considered approxi- 
mately as an LDA exchange-correlation energy for the 
correlated electrons. Because I ex is usually identified as 
the average exchange integral J [TJ, therefore, we con- 
veniently replace the former by the latter in Eq. [13] 
The final expression serves as Ej^ DA in our new LDA+U 
scheme: 
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(14) 



3 



Inserting Eqs. [8] and [l4]in Eq. [2] results in 

U — J 



E' 



E< 



(15) 



It is recognized that the Hartree terms in Ehf and Ef^ DA 
precisely cancel out each other in Eq. [15] 

Accordingly, the corresponding potential for the elec- 
tron of the index n miJ is obtained as the first-order deriva- 
tive of E u with respect to n ma : 

Vma = ~{U - J)n ma . (16) 



When compared with v^ a in Eq. [4] the new poten- 
tial does not contain the extra term, \(U — J), which 
uniformly raises the one-electron potential for all the 
d(f) orbitals. Figure 2 illustrates the spin-orbital en- 
ergy changes when in Eq. 16 is used for a hydrogen 



in Eq. 

atom. In comparison to the one-electron potential in Fig- 
ure 1, the self-interaction correction in v^a- ls complete 
in Figure 2 and the on-site repulsion energy is expressed 
correctly at the same time. 
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FIG. 2: Is spin-orbital energy changes of a hydrogen atom 
by subsequent inclusions of vh = U, {v X c] ~ — J, v xc y ~ 0} 
and {vYsi =-([/- J), vY sl = 0} in the new LDA+U. The 
dashed levels are the unoccupied |-spin energy level after the 
exchange splitting. 

In reality, the total electronic energy is calculated from 
the following expression [5] [8]: 



E total [ Pai h] = E LDA (e 



LDA+U 



)-E u {n). (17) 



-Elua(c^ D ' 4+C/ ) is the total energy component that is ob- 
tained from the LDA expression by using the Kohn-Sham 
orbital energies, e^ Dj4+ ' 7 , calculated from self-consistent 
LDA+U procedures. 

The new LDA+U scheme was implemented in FPLO 
(Full Potential Local Orbital) program (version 5) [8] [9] 
and was tested for the bulk Gd metal. A detailed analy- 
sis of the LDA results on the bulk Gd has been presented 
upon employment of the linearized augmented-plane- 
wave (LAPW) method in a scalar relativistic approxi- 
mation [10]. The unoccupied spin-minority 4/ states are 
centered just 0.5 eV above the Fermi energy, and the 



consequent hybridization of the 4/ states with the de- 
localized 5d states explains the complicated Fermi sur- 
faces observed in experiments and also the appreciable 
contribution of the 4/ electrons at the Fermi energy. 
The low-lying spin-minority 4/ states have been observed 
later in the electronic structures of lanthanide metals and 
their compounds calculated with the SIC-LDA method, 
which accounts for the observed valence states across the 
lanthanide metal series [J]. However, the LDA fails to 
predict the FM ground state of the bulk Gd. It was 
shown that LDA+U (the AL scheme) yields a correct 
FM ground state and this result has been attributed to 
an upward shift of the spin-minority 4/ states from the 
Fermi energy in the LDA+U energy structure [5] ITT] . 
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FIG. 3: the electronic band structures of the FM Gd metal 
from the LDA (top), the LDA+U(AL) (middle) and new 
LDA+U (bottom) methods at the experimental structure. 
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In our LDA+U calculations of the bulk Gd, we em- 
ployed the Perdew-Wang LDA exchange-correlation func- 
tional |12j . For its hep structure, we chose the experi- 
mental lattice constants with a — 6.8662 a.u. and c/a — 
1.587 used in the previous LDA+U(AL) LAPW calcula- 
tions [B]. U = 6.7 and J = 0.7 eV were used for the 4/ 
orbitals (F = 6.70, F 2 = 8.34, F 4 = 5.57, and F 6 = 4.13 
eV) OUT]. For the FM (AFM) state, 133 (372) fc-points 
in the irreducible part of the hexagonal BZ were used, 
with Gaussian smearing. The compression parameters 
were optimized separately for the FM and AFM states 
for 4/, (5s, 5p), 5d, 6s and 6p numerical atomic orbitals. 

Figure 3 shows the FM electronic structures of the bulk 
Gd calculated by LDA, LDA+U (AL) and new LDA+U 
in this work. The LDA result is in an excellent agreement 
with the previous LAPW calculation [10]. The 4/ states 
are separated into two manifolds by an exchange splitting 
of 5.0 eV and thus I ex = 0.71 (= ^) eV [7], conforming 
to the J value employed for the LDA+U calculations [5]. 
In comparison to the LDA, the LDA+U(AL) band struc- 
ture shows a large shift of spin-minority 4/ states by 2.5 
eV, while the spin-majority 4/ states are shifted down- 
wards by 3.5 eV. The effective exchange splitting is larger 
than the LDA result by 6.0 eV, which is consistent with 
the previous LDA+U(AL) LAPW calculations [5]. In 
contrast to the LDA+U(AL), however, the new LDA+U 
scheme does not change the location of the spin-minority 
4/ states, and this leaves the band structure essentially 
unaffected near the Fermi energy. The close proximity of 
those 4/ state around the Fermi energy is consistent with 
the observations from the SIC-LDA energy band struc- 
tures of lanthanide metals and their compounds |3]. The 
effective exchange splitting is again 6.0 eV larger than 
the LDA, which is solely due to the downward shift of 
the spin-majority 4/ states. 

The spin magnetic moment and density of states for 
FM and AFM Gd are listed in Table I for the three cal- 
culations. Both of the spin magnetic moments from LDA 
and the new LDA+U results show a good agreement with 
experiment. The densities of states from the LDA and 
the new LDA+U are much closer to the experiment in 
comparison to the LDA+U(AL) result. Importantly, the 
new LDA+U predicts the correct FM ground state, al- 
though the stabilization against the AFM state is not as 
large as what we obtained from the LDA+U(AL). The 
stability of the FM ground state in the new LDA+U re- 
sults indicates that the upward shift of the spin-minority 
4/ states may not be the origin of the stable FM state, 
in contrast to the previous suggestion (51 ITT]. 

In summary, we have proposed a new LDA+U scheme 
that unerringly corrects the self interaction in the LDA. 
It is based on an MF expression of the sum of the ex- 
change integrals and its connection to the LDA spin- 
polarization energy of atoms. The new method provided 
improved ground-state properties of the bulk Gd, while 
placing the upper Hubbard 4/ states in the energy re- 



gion found as in the SIC-LDA results. Systematic stud- 
ies on various systems will reveal in the future how the 
new method fares. Its close relationship to the existing 
LDA+U schemes allows easy implementation in other 
band structure calculation codes as well. By employ- 
ing an LDA exchange-correlation energy of the correlated 
electrons in the Ej^ DA , the new scheme may eventually al- 
low parameter-free LDA+U calculations because the Eh 
term in the Ef^ DA in our scheme is exact and becomes 
cancelled out in the E u . 
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TABLE I: The spin magnetic moments (in /j,b) and density 
of states at the Fermi energy (in (Ry atom) -1 ) for FM and 
AFM states of the hep Gd from different calculation methods 
(AL: LDA+U(AL); SIC: the new LDA+U in this work). 





LDA 


AL 


SIC 


Expt. 


FM 


7.610 


7.739 


7.657 


7.63 a 




24.17 


14.75 


19.23 


21.35 b 


AFM 


7.384 


7.628 


7.439 






22.09 


17.85 


21.20 




AE C 


-18.9 


37.1 


9.5 





a Reference [13]. 
b Reference 

c S(AFM) - E(FM) in meV/atom. 
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